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ABSTRACT 

More reliable constraints on the microlensing optical depth comes from a 
better understanding of the Galactic model. Based on well-constrained Galac- 
tic bulge and disk models constructed from survey observations, such as, HST, 
2MASS, and SDSS, we calculate the microlensing optical depths toward the 
Galactic bulge fields, and compare them with recent results of microlensing sur- 
veys. We test x 2 statistics of microlensing optical depths expected from those 
models, as well as previously proposed models, using two types of data: opti- 
cal depth map in (/, b) and averaged optical depth over the Galactic longitude 
I as a function of the latitude b. From this analysis, we find that the Galactic 
bulge models of 2MASS, Han & Gould (2003), and G2 of Stanek et al. (1997) 
show a good agreement with the microlensing optical depth profiles for all the 
microlensing observations, compared with E2 of Stanek et al. (1997). We find, 
on the other hand, that models involving an SDSS disk model produce relatively 
higher \ 2 values. It should be noted that modeled microlensing optical depths 
diverge in the low Galactic latitude, |fe| < 2°. Therefore, we suggest the mi- 
crolensing observation toward much closer to central regions of the Galaxy to 
further test the proposed Galactic models, if it is more technically feasible than 
waiting for large data set of microlensing events. 

Subject headings: gravitational lensing — Galaxy: bulge — Galaxy: stellar con- 
tents — Galaxy: structure — methods: statistical 
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INTRODUCTION 



Microlensing survey was originally proposed as a too l for detecting m assive astronomical 

986). Searches for the 



compact halo objects (MACHOs) in the Galactic halo (jPaczyhski 



microlensing e yents towards the Mag ellanic Clouds by the MACHO ( Alcock et al.lll993l ) and 
EROS groups (lAubourg et al.lll993l ) have placed constraints on the fraction of the Galactic 
MACHO populations. Searches for dark objects in the halo of M31 have also been per- 
formed (AGAPE, Ansari et al. 1999; WeCAPP, Riffeser et al. 2003; MEGA, de Jong et al. 
2004; VATT-Columbia, Uglesich et al. 2004; POINT-AGAPE, Calchi Novati et al. 2005; 
Angstrom, Kerins et al. 2006). In addition, microlensing has proven its elf as a powerful tool 



in constraining the distribution of faint stellar object s in the Milky Way (IKiraga fc Paczvriski 



1994; Han fc Gould Il995 



2003 



Wood fe Maoll2005l : IWoodll2007l : ICalchi Novati et al.ll2008l ). 



For instance, lHan &: Gouldl (ll995l ) examined theoretical Galactic models by constructing 
a microlensing map. They explored a triaxial bar-shaped bulge model, w hich was based 



on b y COBE-DIKBE multiwavelength observations of the Galactic bulge ( IWeiland et al. 
19941 ). and demonstrated that observationally determining the microlensing optical depth 
may provide supplementary information on the Galactic mass distribution. 

Unfortunately, however, reality was not so simple. Early measurements of the mi- 
crolensing optical depth towards the Galactic bulge were significantly higher than predic- 
tions, suggesting both observational and theoretical sides make ef forts to reconcile. Mea- 
surements at Baade's window by OGLE from 9 microlensing events ( jUdalski et al.lll994T) and 
at (7,6 ) = (2?55, — 3?64) by MACHO from 13 events of clump giant sources ( Alcock et al. 
1997ah have yielded r = 3.3 ± 1.2 x 10" 6 and r = 3.9±i;l x 10~ 6 , respectively, klcock et al" 



1997a|) pointed out a possibility of a systematic bias such as blending effect in the mi- 
crolensing optical depth measurement to explain those high opt ical depth values in the sense 
that bulge fields are high dense regions. iPopowski et al.l (120011 ) suggested that the bias due 
to the blending can be avoided by using the events only with bright source stars such as 
red clump giant stars. Such stars are identified by their well-defined position in the color- 
magnitude diagram. This position ensures that they are most likely in the Galactic Bulge. 
Their large flux also makes blending problems relatively unimportant. Measurements of 
t using red clump giants have ten ded to reduce disc repancies with models. Based on 16 
events, the EROS-2 collaboration (lAfonso et al.ll2003l ) gave the microlensing optical depth 
of 0.94 ± 0.29 x 10~ 6 at (I, b) = (2?5, -4?0). The MACHO group calculated t he microlensing 



optic al depth toward the Galactic bulge using their seven-year survey data (IPopowski et al. 
20051 ). The measurement of r = 2.17^ x 10~ 6 at (I, b) = (1?50, -2?68) was obtained from 



42 microlensing events detected during the monitoring of about 7.4 x 10 5 clump giant sources 
covering 4.5 deg 2 . More recently, OGLE-II presented the measurement of the microlensing 



optical depth toward the Galactic bulge based on recent four-year survey data (ISumi et al. 
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20061 ). Using the sample of 32 microlensing events in 20 bulge fields covering ~ 5 cleg 2 , 
they found r = 2.55±°;^ x 10~ 6 at (/, b) = (1?16, -2?75). Efforts in a theoretical side also 
refined the microlensing optical depth estimate by using more sophisticated models based 
on observational results from wide field surveys. Theoretical estimates involving a bar in 
the bulge oriented along the line of sight approach to a value in the range of 0.8 x 10~ 6 to 
2.0 x 1Q- 6 (see e.g. IPaczvriski et al.l ll994l; Izhao et al.lll995l : Izhao k Maolll996l : iBinnev et al 



2000l ; lHan & Gouldl 120031 ) . 



The optical depth in microlensing is defined as the probability that a source star is being 
placed within the Einstein radius of a foreground lens star. If bulge stars are distributed 
over a distance d from the Earth, the microlensing optical depth is given by 



4ttG Jo \lo Ds dD i P(A) D] dD s Dl n(D 



r 



J*dD s D*n(D s 



where n(D s ) is the number density of source stars, p(A) is the mass density of lenses along 
the line of sight, Di, D s and Di s are the distances from the observer to the lens and source, 
and the distance from the lens to the source, respectively, and D is DiDi s /D s . 

Given that the microlensing optical depth is subject to an underlying mass distribution 
model, it is crucial to employ a mass model, which can be constructed from independent 
observations, e.g., star counts. In fact, the stellar contents of the disk and the bulge have 
been measured by survey projects and modeled to estimate the optical depth. For in- 
stance, star counts from Hubble Space Telescope (HST) provided constrain ts on faint stars 



down to the hydrogen burning limit in the disk and 0.15 M in the bulge (IHoltzman et al. 



1998; 


Zoccali et al. 


2000; 


Zheng et al. 


2001) 



relatively stron ger constraints on t he mass distribution in terms of the microlensing opti- 
cal depth (e.g 



Han fc Gould 2003 



, , , Results from large-scale sky surveys such as 2MASS 

( Skrutskie et al.lll997 ) and SDSS (jYork et al. 2000 ) are also available in the near-infrared 



and optical wavelength bands. On th e basis of those results, new models of the Ga lactic 
bulge and disk have been suggested (ILopez-Corredoira et al.l 120051 ; lJuric et al.l 120081 ). In 
the present paper, we calculate the microlensing optical depth with Galactic bulge and disk 
models recently by various models based on recent survey observations as well as with those 
studied earlier. We then compare the optical depths based on models with the results of 
MACHO, OGLE-II, and EROS-2 surveys. 

The paper is organized as follows. We begin with a brief description of models of the 
Galactic bulge and the disk in § 2. We present expected microlensing optical depths from 
various models and compare those with observed optical depths in § 3. Finally, we conclude 
with a summary of our results and discussion. 
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GALACTIC BULGE AND DISK MODELS 



2.1. BULGE MODEL 



We calculate theor etical optical depths bas e d on 5 diffe rent mass distribution models. 
To begin with, we follow lLopez-Corredoira et al. J2004! . l2005h . who have suggested two bulge 
models, i.e., triaxial and boxy models, on the basis of the number density from 2MASS star 
counts (ISkrutskie et al.lll997l ). That is, our first and second models are 



PtriaxialV^j Vi z ) 



p T exp 



and 



respectively, where 



pboxy(x,y,z) = p B exp 



740 pc 
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(2) 



(3) 



(4) 



(5) 



The normalization in the triaxial and boxy bulge models are px = 9.9 M /pc 3 and pe = 4.67 
M /pc 3 , respectively. Following ICalchi Novati et al.l (120081 ) we assume the total bulge mass 
of 1.5 x 10 10 Mq within 2.5 kpc from the Galactic center. Note that this value is in a range 



of generally accepted bulge ma ss of the Galaxy 1 — 2 x 10 10 M Q (e.g. lBlumlll995l ; IZhao et al. 
19961 ; iDehnen fc Binneyl Il998l ). For both bulge models, the longest axis is inclined by 29° 
with respect to the line of sight. 

We also e mploy three w i dely-u sed bulge models. Our third bulge model (HG03) is 
one studied by Han fc Gouldl (20031). which is favored by an analysis of the COBE DIRBE 



observations (IDwek et al.lll995l ) 



Han &: Gouldl (120031 ) normalized the model u sing the HST 



star c ount. For the fourth model, we adopt G2 model, which can be found in IStanek et al. 
(119971 ). given as 

(6) 



p G2 (x,y,z) 



Poexp 



where 




-| 1/4 



(7) 
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Following IStanek et al.l (119971 ) we take xq = 1239 pc with the axis ratio of xp : yn '■ Zn 
2.8 : 1 .4 : 1 and an inclination angle a = 24?9. The last model is E2 model of IStanek et al. 
(119971 ). The density distribution is given in the form of 



p E2 {x,y,z) = p exp(-r), 



(8) 



where 



x + y 



z 



-I 1/2 



(9) 



Here we take xq = 897 pc with the axis ratio of xq : yo ■ zq = 3.6 : 1.5 : 1 and an inclination 
angle of the bulge major axis with respect to the line of sigh t a = 23?8. For HG03, G2 , and 
E2 bulge mod e ls, we follow the bulge mass normalizations of ICalchi Novati et al.l (120081 ) and 
Han fc Gouldl (120031 ). We set the distance to the Galactic center as 8.0 kpc throughout our 
analysis. 



2.2. DISK MODEL 



For 2MASS bulge models we employ two disk models (dis 



derived from 2MASS star counts (ILopez-Corredoira et al. 



i 1 mo del and disk 2 model) 



2004 



2005|). For bulge models of 



HG03, G2, and E2 we adopt the disk model of IZheng et al.l (120011 ). whose density profile 
is given as a sech 2 (exponential) function for the thin (thick) components. We normalize it 
with the disk mass density of the Solar neighborhood, p = 0.05 M /pc 3 (cf. lHan fc Gould 
20031 ). In addition, we also investigate t he SPSS disk mod el with respect to all the bulge 
models we consider in the present paper. lJuric et al.l (120081 ) estimate the three-dimensional 
number density distribution of the Galactic disk and halo using the photometric parallax 
method on the SDSS data. We use their disk models consisting of thin and thick exponential 
disk with local thick-to-thin disk normalization pthick(-R©)/pthin(-R©) = 12%. Again, the local 
number density is normalized to p(Rq) = 0.05 M /pc 3 . 



3. RESULTS 

In Figures 1 and 2, we show the map of the optical depth difference normalized by the 
observational error a T for each (/, b) field, which is given by e = {r — r m )/a T , where r D and r m 
are the microlensing optical depths in a given field from observation and model, respectively. 
We compare the microlensing optical depth observed from MACHO and OGLE-II with those 
calculated from models in field by field. Figures 1 and 2 result from MACHO and OGLE-II, 
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respectively. In the map, the fields of negative e and positive e are shown in blue and red, 
respectively. For those of r Q = 0.0, fields are marked by black. Dotted curves in each panel 
represent the microlensing optical depth contour of a model in units of 10~ 6 . Each panel 
of Figures 1 and 2 results from different mass models indicated in the lower left corner of 
the panel. Table 1 summarize how to pair a bulge model and a disk model. For instance, 
2MT1 and 2MB1 represent combinations of 2MASS triaxial bulge model and the 2MASS 
disk 1 model, 2MASS boxy bulge model and the 2MASS disk 1 model, respectively. In case 
of 2MT2 and 2MB2, the 2MASS disk 1 model is replaced by the 2MASS disk 2 model for 
a given bulge model. When the SDSS disk model is substituted, it is indicated as shown: 
HG0 3, G2, and E2 st and for cases of HG03, G2, and E2 bulge models and the disk model 
of Izheng et all toOlh . 



To quantify the difference of the microlensing optical depth, we calculate xfby defined 



as 



Xl.b 



N 

E 



Tvt 



(10) 



where r 0j j, r m i and a Qy i are the microlensing optical depths from observation and model, 
and the observational error in each field, respectively. Here summations are repeated over 
14 OGLE-II fields and 19 MACHO fields excluding fields of r G = 0. We summarize results 
in Table 2. As seen, one finds that 2MASS, HG03, and G2 models agree well with the 
microlensing optical depth dist ribution for MACHO and OGLE-II when disk models of 
2MASS and IZheng et al.l (120011 ) are used. It should be noted that models that include 
the SDSS disk model produce relatively higher xfb values. We note that the microlensing 
optical depth from MACHO has much higher xfb va hies than OGLE-II since the microlensing 
optical depth values in four MACHO fields (176, 113,109, and 105) are 3a below with 
respect to the microlensing optical depth values expected from all models in those fields. 
Therefore, it should be pointed out that comparing two columns of OGLE-II and MACHO 
is meaningless. We repeat calculation of \ 2 without those four fields and show results in 
parentheses for comparison in Table 2. Conclusion remains same. The reason why four of the 
MACHO fields deviate out of 3 a could be either that the models are too smooth to explain 
details of Galactic structures while there could be components affecting more optically than 
dynamically or that observational contaminations like the blending effect may play a role. 
For instance, OGLE-II investigates the blending effect even in bright ev ents and points out 



Alcock et al. 


1997b: 


Smith et al. 


2007) 



(e.g. 



Udalski et al. 1994 



In Figure 3, we also show microlensing optical depths averaged over the Galactic lon- 
gitude as a function of the Galactic latitude b. In the left panel, averages of microlensing 
optical depths measured by OGLE-II, MACHO, and EROS-2 are shown with their error bars. 
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OGLE-II, MACHO, and EROS-2 data are plotted with filled circles, squares, and triangles, 
respectively. Best fits are superposed on top of observed data. Different line types represent 
different microlensing experiments: continuous, short dashed and long dashed lines stand 
for best-fit curves of OGLE-II, MACHO, and EROS-2, respectively. Analytic functions of 
those fits estimated by OGLE-I I, MACHO and EROS-2 are given as t OG le-ii = [(4.48 ± 
2.37) + (0.78 ±0 . 84) x b] x f (T 6 JSumi et alil2~0~0~6L t M acho = [5.2 + (f .06 ± 0.71) x b] x 10~ 6 
f Popowski et~aD l2005h. and r ERO s-2 = (1-62 ± 0.23) exp[-(0.43 ± 0.16)(|6| - 3°)] x 10~ 6 
f Hamadache et al. 2006I ). respectively. In the right panel, microlensing optical depths to- 
ward the Galactic bulge expected from models described in the previous section are shown 
as a function of the Galactic latitude b with thin curves. Dashed curve represents the cal- 
culated microlensing optical depth resulted from the 2MASS triaxial bulge model and the 
2MASS disk 2 model (2MT2). Dotted curve represents the calculated microlensing optical 
depth resulted from the 2MASS boxy bulge model and the 2 MASS disk 2 model (2MB2). 
Long dashed, dot-long-dashed, and dot-short-dashed curves represent HG03, G2, and E2 
models, respectively. For comparison, we plot best fits obtained by different microlensing 
experiments with thick curves, which are already shown in the left panel. Continuous, short 
dashed and long dashed curves correspond to results of OGLE-II, MACHO, and EROS-2, 
respectively, as in the left panel. Models shown in the plot in general seem compatible with 
all the microlensing observations within the observational uncertainties. It is interesting to 
note that deviations among model fits become larger as the Galactic latitude decreases. To 
evaluate more quantitatively we calculate xl i n eacn & using equation (10), and list results 
in Table 3. As seen previously in Table 2, 2MASS, HG03, and G2 result in smaller xl m au 
the microlensing experiments. In Table 3 we show xl of all cases, including those not shown 
in the right panel of Figure 3. One may find again that when SDSS disk model is used fits 
become poorer. 



4. SUMMARY AND DISCUSSIONS 

By comparing the microlensing optical depth map constructed by a mass model with 
the observed microlensing optical depth one may put a tight constraint on the Galactic 
model. We have estimated the microlensing optical depths toward the Galactic bulge using 
the Galactic bulge and disk models constructed from survey observations, and compared 
them with the sample of microlensing events observed towards the Galactic bulge with red 
clump giant sources reported by the MACHO, OGLE-II and EROS-2 collaborations. From 
the analysis we carried out, we show that model estimates are in general compatible with 
microlensing observations. According to calculated x 2 we find that 2MASS, HG03, and G2 
models reproduce the microlensing optical depth distribution from observations well. We 
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also find that for those well-fit models contribution to the microlensing optical depth of 
the disk components is not crucial. For example, for 2MASS cases two disk models do not 
make any difference in \ 2 ■ It is, however, interesting to note that the analyses of SDSS 
disk models give us somewhat different results. That is, models including SDSS disk model 
produce relatively higher x 2 values than others do. 

It should be pointed out that current observational data of both the microlensing optical 
depth in two dimensional fields and the averaged optical depth over I are yet statistically 
insufficient to rule out conclusively any specific models we have considered since the number 
of observed microlensing events from monitoring RCGs in recent microlensing surveys is still 
small. On the other hands, we note that the difference of the average optical depth among 
models increases as closer to the Galactic center. Therefore, we suggest the microlensing 
observation toward closer to the Galactic center regions if possible may be more effective 
in constraining the Galactic model than extending microlensing searches to detect more 
microlensing events. 
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This work is the result of research activities at the Astrophysical Research Center for the 
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gineering Foundation (KOSEF). HYC was supported by the Korea Research Foundation 
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Table 1. Bulge and disk models 



2MASS disk 1 2MASS disk 2 Zheng disk SDSS disk 



2MASS triaxial bulge 
2MASS boxy bulge 
Han & Gould (2003) 
G2 of Stanek et al. (1997) 
E2 of Stanek et al. (1997) 



2MT1 2MT2 
2MB1 2MB2 



2MT+SDSS 
2MB+SDSS 
HG03 HG03+SDSS 
G2 G2+SDSS 
E2 E2+SDSS 



Table 2. The chi-square values of microlensing optical depths between observations and 
model estimates in two dimensional microlensing optical depth map. 





Model OGLE-II MACHO 



2MT1 


13 


.73 


122. 


23 


(27. 


20) 


2MT2 


13 


.57 


122. 


31 


(27. 


51) 


2MT+SDSS 


19 


.23 


178. 


20 


(38. 


76) 


2MB1 


15 


.24 


143. 


.84 


(34. 


55) 


2MB2 


14. 


.91 


142. 


.55 


(34. 


63) 


2MB+SDSS 


21 


.84 


208. 


.89 


(49. 


50) 


HG03 


14. 


.90 


121. 


.72 


(24. 


92) 


HG03+SDSS 


17 


.68 


156. 


.83 


(32. 


77) 


G2 


16 


.81 


159. 


30 


(35. 


15) 


G2+SDSS 


20 


,13 


202. 


.00 


(46. 


15) 


E2 


23 


.66 


200. 


.51 


(41. 


89) 


E2+SDSS 


28. 


.05 


248. 


.51 


(54. 


15) 



Note. - The number of fields used in 
the summation is 14 and 19 for OGLE-II 
and MACHO, respectively. For MACHO, 
chi-square values excluding 4 fields whose 
errors are greater than 3 a are also shown 
in parentheses for comparison. 
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Table 3. The chi-square values of averaged optical depths as a funciton of b. 



xl 



Model 


OGLE-II 


MACHO 


EROS-2 


2MT1 


2.27 


1.66 


10. 


.71 


2MT2 


2.26 


1.76 


10. 


.57 


2MT+SDSS 


3.48 


6.94 


22, 


03 


2MB1 


2.76 


4.41 


18. 


.41 


2MB2 


2.73 


4.35 


18. 


.00 


2MB+SDSS 


4.81 


13.41 


35. 


.59 


HG03 


2.41 


1.76 


10. 


.83 


HG03+SDSS 


3.13 


5.20 


17. 


.46 


G2 


3.30 


6.25 


21. 


11 


G2+SDSS 


4.54 


12.05 


31. 


.55 


E2 


4.94 


13.42 


28. 


00 


E2+SDSS 


6.60 


21.26 


39. 


.22 



13 
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Fig. 1. — Relative difference of microlensing optical depths in 10~ 6 between MACHO results 
and models, normalized by the observational error a T , as a function of (l,b). Note that 
the fields of negative e and positive e are shown in blue and red, respectively. For those of 
t = 0.0, fields are marked by black. Dotted curves in each panel represent the microlensing 
optical depth contour of a model in units of 10~ 6 . Each panel results from different mass 
models indicated in the lower left corner of the panel. See the text for the definition of the 
error-normalized difference and for model descriptions. 
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Fig. 2. — Same plots as Figure 1, but showing results from OGLE-II. 
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Fig. 3. — Averaged microlensing optical depths overthe Galactic longitude I as a function of 
the Galactic latitude b. Left: Averages of microlensing optical depths measured by OGLE-II, 
MACHO, and EROS-2 are shown with their error bars. OGLE-II, MACHO, and EROS-2 
data are plotted with filled circles, squares, and triangles, respectively. Continuous, short 
dashed, and long dashed lines stand for the best-fit curves of OGLE-II, MACHO, and EROS- 
2, respectively. Right: Microlensing optical depths toward the Galactic bulge expected from 
models are shown as a function of the Galactic latitude b with thin curves. Dashed and 
dotted curves represent the calculated microlensing optical depth resulted from 2MT2 and 
2MB2, respectively. Long-dashed, dot and long-dashed, and dot and short-dashed curves 
represent HG03, G2, and E2, respectively. For comparison, OGLE-II, MACHO, and EROS-2 
microlensing experiment results are shown with thick continuous, short-dashed, and long- 
dashed curves. 



